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ABSTRACT 

Using mid-IR and optical data, we deduce the total infrared (IR) luminosities 
of galaxies in the Coma cluster and present their infrared luminosity function 
(LF). The shape of the overall Coma IR LF does not show signihcant differences 
from the IR LFs of the general held, which indicates the general independence of 
global galaxy star formation on environment up to densities ~ 40 times greater 
than in the held (we cannot test such independence above Lir ~ 10^^ ergs s“^). 
However, a shallower faint end slope and a smaller L*. are found in the core 
region (where the densities are still higher) compared to the outskirt region of 
the cluster, and most of the brightest IR galaxies are found outside of the core 
region. The IR LF in the NGC 4839 group region does not show any unique 
characteristics. By integrating the IR LF, we hnd a total star formation rate in 
the cluster of about 97.0 We also studied the contributions of early- and 

late-type galaxies to the IR LF. The late-type galaxies dominate the bright end 
of the LF, and the early-type galaxies, although only making up a small portion 
(~ 15%) of the total IR emission of the cluster, contribute greatly to the number 
counts of the LF at Lir < 10"^^ ergs s“^. 

Subject headings: galaxies: clusters: individual (Coma) — galaxies: luminosity 
function — infrared: galaxies 


1. Introduction 

Galaxy evolution is largely the story of how the masses, morphologies, and patterns of 
star formation in these objects vary with environment and cosmological epoch. Luminosity 
functions (LFs) are very important statistical tools for studying evolutionary changes in 
galaxy populations and provide key observational constraints on galaxy evolution. The 
infrared (IR) LF, in particular, by tracing the IR emission of dust heated by star forming 
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activities, remains one of the best probes to study the global evolution of the star formation 
rate (SFR) with environment and redshift. All the indicators show that the SFR is a strong 
function of epoch. The average SFR has declined by an order of magnitude to the present 
epoch from a peak near z 1 (Lilly et ah 1996; Madau et ah 1996; Perez-Gonzalez et ah 
2005; Le Floc’h et al. 2005). 

We also expect the SFR to depend on environment. In crowded regions, e.g., the 
cores of dense clusters, galaxies experience many interactions and mergers, which drive the 
accumulation of mass and redistribute the gas. The hierarchical galaxy formation models 
(Somerville & Primack 1999; Cole et al. 2000) suggest that galaxies in a cluster gradually 
lose their gas reservoir as they are accreted into the cluster center. Those environmental 
effects are expected to result in a systematic modihcation of the patterns of star formation 
in clusters. In fact, a number of studies (Gomez et al. 2003; Balogh et al. 1998) have shown 
such effects, with a lower SFR in clusters compared with the held. The infrared LFs of 
clusters, as a global indicator of SFR, may also carry such an imprint and differ from the IR 
LFs of held galaxies. 

Since the launch of the Infrared Astronomical Satellite (IRAS) and the Infrared Spaee 
Observatory (ISO), there have been many studies of the infrared LF of galaxies (e.g.. Rush, 
Malkan & Spinoglio 1993; Serjeant et al. 2001; Takeuchi, Yoshikawa & Ishii 2003). However, 
these works concentrate on general held surveys and there is no systematic study of IR LFs 
with epochs devoted to the cluster population. Such a study will enable us to disentangle 
the evolution of the SFR in diherent environments and help us understand where and why 
the changes in SFR occurred. 

The recently launched Spitzer Space Telescope with the Multiband Imaging Photometer 
(MIPS, Rieke et al. 2004), capable of high resolution, large sky surveys with high sensitivity, 
provides us the capabilities necessary for thorough studies of the star forming properties of 
dense galaxy clusters at different redshifts. In this paper, the hrst one of a series studying 
the IR LF of clusters up to 2 : ~ 0.8 using MIPS 24 and 70 /rm observations, we present 
the IR LF of the Goma cluster. The Goma cluster, as the nearest rich cluster, provides 
us an excellent chance to study the properties of the IR LF down to a very faint limit. It 
also enables us to study the change of IR LF in different regions of the cluster, as well as 
the contributions to the total IR LF from different types of galaxies. This paper will show 
how the cluster environment shapes the current star formation in this prototypical dense 
cluster and it will provide a foundation for the future studies of the SFR patterns in other 
clusters, both nearby and at high redshift. In the paper, we use the cosmological parameter 
set (h, Go, Ao) = (0.7, 0.3, 0.7). We assume a distance modulus of m - M = 35.0 mag for the 
Goma cluster at z=0.023 (Struble & Rood 1999). 



3 


2. Observations and Data Reduction 

We used MIPS to observe the Coma cluster in medium scan map mode on Jun 6, 2004. 
Our map covered a 4 deg^ area centered at a = 12^59’”27^, S = +27°51'53", which included 
both the cluster core region and the NGC 4839 region. Fig. 1 shows the region of the 
map. The 24 /im and 70 /im reductions were carried out with the MIPS Data Analysis Tool 
(Gordon et al. 2005). The total exposure time was about 88 seconds per pixel for the 24 yum 
observations and about 40 seconds for the 70 /im ones. 


3. Source Detection 

SExtractor was applied to the images to detect sources automatically and to obtain 
photometry parameters. First, the images were analyzed and sky background models were 
built. Then the images were background-subtracted and hltered with Gaussian functions 
with the full width at half maximum (FWHM) matching the FWHM of the MIPS 24 /im 
and 70 /im point spread functions (PSFs). All the objects with values exceeding a certain 
threshold of local background were extracted. We set this detection threshold, relative to 
the root-mean-square background noise, at 0.65 for the 24 /im image and 2 for the 70 /im 
image. After deblending the adjacent objects and ‘cleaning’ the artifacts due to bright 
objects, SExtractor gave a hnal catalog of sources. SExtractor provides several types of 
magnitude measurements. We adopted the MAG_BEST magnitude. In most of the cases, 
the MAG_BEST magnitude is measured in an adaptive aperture of 2.5 times the Kron 
aperture, but in crowded helds, it is measured in an isophotal aperture and corrected for 
aperture losses. In our data, few regions are crowded. 


4. Completeness 

The completeness of source detection affects the LE directly, so it is important to know 
the detection limit of the observations. Papovich et al. (2004) studied the source detection 
completeness of MIPS 24 //m images in several different helds. Among these helds, the 
Bootes image has almost the same exposure time as the Goma image and the background 
levels are also similar, with a mean value of 22.7 MJy sr“^ for the Bootes held and 33.4 MJy 
sf-i fQr Goma held. By inserting artihcial sources in the images and performing source 
extraction on them, Papovich et al. (2004) found a 80% completeness hux density limit of 
0.27 mJy in the Bootes held. As a rough approximation, a simple linear scaling with the 
square root of the background level gives us a detection limit of 80% completeness at 0.33 
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mJy for the Coma field. 

At 70 /im, the Coma data were obtained at a lower detector bias than those in Bootes, 
resulting in an improvement in overall performance. Therefore, we ignore the background 
difference and adopt the completeness limit of about 80 mJy obtained in the Bootes held 
with the same exposure time (Dole et ah 2003). 


5. Matching Spectroscopic Surveys of the Coma Cluster w^ith 24 fim. Sources 

To study the infrared LF of the Coma cluster galaxies, we used a spectroscopic sample so 
that cluster membership could be conhrmed. Fortunately many spectroscopic surveys have 
been carried out in this region. Among them, the catalogs from Beijersbergen & van der 
Hulst (2003) and Mobasher et ah (2001) have the largest overlapping area with the MIPS 24 
/im observations and also go to fairly deep detection limits. Beijersbergen & van der Hulst 
(2003) generated a catalog (hereafter BvdHC) using all known Coma cluster redshifts in a 5.2 
deg^ region. This catalog covers almost the whole region of the MIPS 24 /rm image except 
a few small patches at the edge and has a 93% completeness down to Sloan r' = 16.27 mag. 
Mobasher et al. (2001) performed spectroscopic observations on two rectangular regions, one 
at the cluster core and the other near the NGC 4839 group, each 32.5 x 50.8 arcmin^. The 
core region is totally covered by MIPS observations, and the region near the NGC 4839 
group is partly covered. Its completeness is about 60% for the bright galaxies (i? < 17) and 
decreases towards the faint end. The difference between the Sloan r' hlter BvdHC used and 
the Cousins R hlter Mobasher et al. (2001) used is small. Comparing the common objects in 
these two catalogs gives a difference of R — r' ~ 0.03 mag, so we do not differentiate them and 
just use R to refer to both of them in this paper. Furthermore, we complemented Mobasher’s 
catalog with the BvdHC down to R = 16.27 mag, and generated a merged catalog in these 
regions (hereafter MBC). 

We selected all the galaxies from these two catalogs (BvdHC & MBC) with 4000 
km s“^ < cx < 10000 km s~^ as cluster members (Colless & Dunn 1996) and cross matched 
them with our 24 /xm and 70 /rm sources. Any 24 /rm or 70 /rm source within 10" of the 
optical galaxy was identihed as the IR emission from this galaxy. This search radius is about 
twice as large as the FWHM of the MIPS PSF at 24 /im and half of the FWHM at 70 
yum. It allows a displacement in projected distance up to 5 kpc between the optical centers 
of galaxies and the peaks of their IR emission. About 90% of the cluster members with 
24 yum emission above the completeness limit have a displacement between optical center 
and IR emission peak smaller than 5", i.e., 2.5 kpc. When multiple identihcations occurred, 
the one with the smallest distance from the optical galaxy was selected. Less than 2% of 
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the sources had multiple identihcations. Therefore, our dual sample is not sensitive to the 
details of matching infrared and optical sources; few cases yield ambiguous associations, and 
moderate changes in the acceptance radius have little effect on the results. 

Among the 498 Coma galaxies in the BvdHC within the MIPS 24 /im field, 217 have 
24 /im counterparts. In the part of the held covered at 70 /rm , 58 were detected out of 477 
members. In the MBC, there were 123 galaxies detected at 24 /rm out of 333 galaxies and 33 
at 70 //m out of 302 galaxies. The number of galaxies detected in both bands was 56 for the 
BvdHC and 33 for the MBC. Although the total number of galaxies detected in the MBC 
is less than in the BvdHC, the overlapping area with MIPS observations is also smaller; it 
is about 0.8 deg^ for the MBC and about 3 deg^ for the BvdHC. Therefore, the number 
density of galaxies detected in the MBC is still larger than in the BvdHC, consistent with 
their different detection limits. 


6. IR Luminosity Function 
6.1. Determination of total IR luminosity 

Since we detected relatively fewer galaxies at 70 /rm than at 24 /rm, we based our LF 
calculations mainly on 24 /rm sources. 

To obtain the total IR luminosities of galaxies, which relate to the total flux from 8 - 1000 
/im, a single measurement of flux density at 24 /rm is not enough. We need more constraints. 
Using a self-consistent modelling of the spectral energy distributions (SEDs) of galaxies over 
a broad range of wavelength, Devriendt, Guiderdoni & Sadat (1999) published a sequence 
of galaxy SEDs with different IR luminosities based on a sample of nearby galaxies. Their 
sample includes normal spirals, luminous IR galaxies (LIRGs) and ultraluminous IR galaxies 
(ULIRGs). If we assume these SEDs are a complete assembly of representatives of nearby 
galaxies, the color correlation of these SEDs should be the same as the color correlation of 
the Coma galaxies. More specihcally, if we know the color correlation between the ratio of 
IR Inminosity (Tjr) and 24 /rm Inminosity (T 24 ) and the ratio of flux density at 24 /xm (S'24) 
and R band {Sn) from the template SEDs, we will know the color correlation of the Coma 
galaxies as well. Therefore, we can use observational data regarding L 24 , Sr and S '24 to get 
the total infrared Inminosity Lr. 

However, the assnmption that the template SEDs include all the galaxy types in Coma 
is not correct. For a cluster as rich as Coma, the early type (E/SO) galaxies dominate 
the optical emission of the cluster. In the infrared, the spiral galaxies are generally more 
Inminous than the early type galaxies and hence are the majority of the IR sonrces. However, 
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given the sensitivity of MIPS and the closeness of the Coma cluster, we still detected the 
24 /xm emission from many elliptical galaxies and SO galaxies. In the BvdHC, which gives 
information on the galaxy type, about half of the galaxies detected at 24 pm are early type 
galaxies, and the rest are mostly spiral galaxies or galaxies without type identihcation. The 
infrared emissions of the early-type galaxies may come from different physical mechanisms 
or different dust geometry than that of the spirals, and their SEDs may have different shapes 
and colors than those of the template SEDs. 

To check for possible differences between the early-type galaxies and the spiral galaxies 
in the colors that are crucial to determine the total infrared luminosity, we plot in Fig. 2 
the color correlation of the different types of galaxies in the BvdHC. Panel a in Fig. 2 is 
the color-color plot of S' 7 o/S '24 vs. S 24 :/Sr for the galaxies detected both at 24 pm and 70 
pm and with a morphology identihcation. The plot shows that early type galaxies (open 
circles) have smaller 824 ,/Sr ratios on average than the spiral galaxies (open triangles), but 
their S' 7 o/S '24 ratios are similar to the spiral galaxies. Although with a large dispersion, the 
template SED’s (crosses) color correlation represents the average value for the whole galaxy 
sample fairly well. If the IR emission of the galaxy mainly comes from dust at a single 
temperature, then small differences in the 810/824 color indicate a small difference in the 
Lir/L 24 ratio for these galaxies. In panel b, we show the 24 pm hux density vs. 824 /Sr 
color of all the galaxies. The early type galaxies mostly reside in the lower corner of the plot 
but they are well mixed with the spiral galaxies and show no difference in this correlation 
from the faint spiral galaxies. The plot also shows that most of the galaxies with 24 pm flux 
density larger than 6 mJy are also detected at 70 pm (indicated by the hlled symbols). This 
result is consistent with the detection completeness at 70 pm assuming the average S 10 /S 24 
color. There also appears to be a trend between the 24 pm flux and the S 24 /SR color of 
the galaxies, which indicates larger IR emission from the redder galaxies. Using this trend, 
we can also obtain a detection limit set by the completeness of the 70 pm observation for 
panel a shown as the dotted line; the region at the left and below the line is affected by the 
incompleteness and large uncertainties of the 70 pm measurements. From both plots, we 
hnd that the early type galaxies are generally less luminous at 24 pm than the spiral galaxies 
and therefore have lower 824/877 ratios, but their color correlation does not differ from that 
of the spiral galaxies with similar 24 pm flux density. This justihes the use of the template 
8 EDs as a complete assembly of all types of galaxies to deduce the total IR luminosity from 
the Lir/L 24 vs. 824 / Sr correlation. 

To obtain the flux densities of the 8 EDs at different bands, we convolved the 8 EDs with 
the response functions of the hlters. For the 24 pm and 70 pm bands, we also account for the 
color corrections as described in the MIP 8 Data Handbook. In Fig. 3, we plot the correlation 
of Lir/L24 vs. 824/Sr obtained from template 8 EDs as well as the value interpolated from 
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the Coma galaxies. The log(5*24/S'i?) ratios of the template spirals range from about —0.5 
to 1 and those of the LIRGs and ULIRGs from about 2 to 3.5. Most of the Coma galaxies 
have log(5*24 /*S'r) ratios smaller than 1.5 and only one source has a color similar to the 
LIRGs/ULIRGs. The \og{Lij./L 2 i) ratios of the LIRGs/ULIRGs are almost a constant of 
~ 1.5. These ratios for the normal spirals increase slowly with the decrease of 824 /Sr 
ratios, with a little (and insignificant) dip at about log(5*24/S'j?) ~ 0. Although a simple 
interpolation onto the correlation works for many of the Coma galaxies, our template SEDs 
do not cover the range with log(5'24/>S'_R) < —0.5 as the data do, so the color correlation at 
this end is an extrapolation from the last few points of the template SEDs. This unbounded 
extrapolation may cause systematic errors when deducing the from the Ljr/L 24 ratio. 
However, we can End some support for the higher value of the 8^/824 ratio at this end from 
panel a of Eig. 2. Despite the incompleteness and large uncertainties, the panel shows a 
slightly higher value of the S 10 /S 24 ratio than the ratio given by template SEDs towards the 
lower end of 82 ^/ Sr ratio. It is worth noticing that this extrapolation is also consistent with 
the general expectations for thermal emission: galaxies with smaller values of 824 / Sr have 
relatively less emission by warm dust and therefore a higher value of Li^j 824 ,. In the future, 
a template SED with lower 82 ^/ Sr ratios will be needed to further constrain the Ljr/L 24 
color at this end, but for now we rely on the simple extrapolation. With the 8^1824 ratio 
of each galaxy in hand, we can directly deduce the 8 ir of each galaxy from its 824 . 

An important result from Fig. 3 is that 824 ,/ Sr is nearly independent of 8 irl 824 over 
the luminosity range of interest to us. Therefore, our initial galaxy selection on the basis of a 
visible spectroscopic study will not introduce biases in the infrared properties of the sample. 

To test the method we used here to determine 8 ir, we compare our result with the 
IR luminosities obtained with the method used by Le Floc’h et ah (2005). They adopt a 
different set of SED templates that are luminosity dependent, e.g., the templates given by 
Lagache et ah (2003). Lagache et al.’s galaxy templates include separate SEDs for normal 
galaxies and starburst galaxies. Their normal galaxies, again, only include spiral galaxies. 
As we can see from Fig. 3, most of the galaxies in the Goma cluster are normal galaxies and 
only one galaxy has a Sr /824 color similar to LIRGs/ULIRGs, so we plot Lagache et al.’s 
correlation between 81^/824 for normal galaxies in Fig. 3 as the dashed line. This correlation 
agrees well with the 8 ir /824 ratios given by the template spirals and thus demonstrates the 
consistency between our method and that of Le Floc’h. 

To test the self-consistency of our method, we use the data for the members of the 
BvdHG detected both in the 24 /rm and 70 pm bands. Using the same SED mapping 
method as before, but, using different combinations of bands, we estimated new values of 
8 ir for these objects. We plot them against the previous obtained from the correlation 



between Ljr/-b 24 and S 24 /SR in Fig. 4. The top panel in Fig. 4 is the Lir obtained from 
the correlation between and S^iq/Sr compared with the previous one. The bottom 

panel is the obtained using the correlation between Lir/LjQ and S-ro/S 24 . The Lir values 
obtained from different color correlations are generally consistent with a standard deviation 
of about 0.10 dex for the top panel and 0.16 dex for the bottom panel. The galaxies with 
S '70 under the 80 mJy completeness limit, shown as the open circles in the figure, have large 
uncertainties in their Sjq measurement (with errors up to 40 %) and therefore show a more 
scattered correlation. The dispersions of the correlations are generally consistent with the 
dispersion caused by the uncertainties in the photometric measurements. It is also possible 
there are signihcant contributions to the dispersions by intrinsic color dispersions of the 
galaxies, as opposed to the tight correlation we assumed in Fig. 3. Since the dispersions are 
modest and have zero averages, they will in any case have little effect on the LF we deduce. 


6.2. Contamination from AGNs 

When we measure the IR emission from galaxies and study their star forming activities, 
contamination from AGNs is always an issue. The IR emission of AGNs comes from dust 
heated by the soft X-ray and ultraviolet emission of the active nuclei rather than from star 
forming activities; therefore, their SEDs could be very different from the template SEDs we 
used. 

To search for the AGNs in the Goma cluster, we cross matched the Quasars and Active 
Galactic Nuclei catalog (Veron-Getty & Veron 2003) with the BvdHG and MBG. The Veron- 
Getty & Veron (2003) catalog is not complete but it includes almost all the AGNs in the 
literature. There are three AGNs detected at 24 /rm : D 16, Syl; NGG 4853, Sy and KUG 
1259-1-280 Sy. Among them, only NGG 4853 has a Lir greater 10^^ ergs s“h These AGNs, 
so few in number, do not have a noticeable effect on the IR LFs we obtain. 


6.3. Total IR Luminosity Function 

After testing for the method we used to deduce the total IR luminosity as described 
above, we calculated the projected IR luminosity function of the Goma cluster. 

For the BvdHG, we obtained the number density of galaxies per projected area by 
directly counting the number of galaxies in each luminosity bin and dividing the number 
by the projection area. For the MBG, we assume the completeness function is unity for 
R < 16.27 mag and behaves as described in Mobasher et al. (2001) for galaxies fainter than 
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R = 16.27 mag. In calculating the number counts of galaxies, we used the inverse of this 
completeness function as a weighting factor to correct for the incompleteness. 

Both LFs are affected by the completeness of the spectroscopic surveys as well as the IR 
observations. The BvdHC spectroscopic completeness is about R = 16.27 mag, and for the 
MBC, since we already correct for the incompleteness, the completeness limit can be extended 
to 7? = 19 mag. The spectroscopic completeness levels of the BvdHC and the MBC are both 
well above the detection limits in surface brightness for the two surveys (Beijersbergen et al. 
2002a; Komiyama et al. 2002), suggesting our samples are not limited by the galaxy surface 
brightness. In the end, we need to estimate the IR completeness of samples selected in the 
optical for the spectroscopic surveys. We therefore need to relate the two spectral ranges. 
About 93% of our galaxy sample has a 824 ,/Sr ratio smaller than 6.5, so we can use this 
value to set the upper limit of the 24 /rm flux density corresponding to the completeness in 
R band. With the linear correlation between Lir and L 24 given by Lagache et al. (2003), 
we hnd the completeness limits of the two spectroscopic surveys in total IR luminosity are 
2.6 X 10® Lq and 2.2 x 10® L©, respectively. From section 4, we already know that the 80% 
detection limit of the 24 /rm observations is about 0.33 mJy, which corresponds to a total 
IR luminosity of about 1.41 x 10® Lq. This detection limit is lower than the completeness 
levels set by spectroscopic surveys. Thus, the spectroscopic data are the primary limit on 
the range at which the LFs are free from the effects of incompleteness, with only ~ 7% 
incompleteness at the lowest luminosities from the dispersion in IR properties. 

In the above calculations of the completeness, we utilize the S24/Sr color distribution 
of our Coma galaxy samples. However, an issue with these samples is that they are optically 
selected and therefore they are likely to miss galaxies faint in the optical while bright in the 
IR. The IR/optical color distribution of our samples might be tilted towards lower values 
and we might underestimate the number of galaxies we missed in the calculation of the com¬ 
pleteness limit. To check the color bias of our optically selected samples, we use the catalog 
of Karachentsev (2004), which provides a nearly complete listing of local galaxies within 10 
Mpc. The catalog gives the R-band magnitude of each galaxy. We obtained the IRAS data 
for these galaxies, and therefore have a complete sample which is not constrained at the 
levels of interest by the optical and/or IR detection limits. We calculated the IR/optical 
colors 824/83 of the normal galaxies (not dominated by AGN and below the LIRG lumi¬ 
nosity range) in the sample and compared them with those of the BvdHG. In B band, the 
BvdHG is complete down to R = 17.5 mag, and the IR completeness limit of 2.6 x 10® Lq 
corresponds to a 824/83 color of ~ 14.6. About 95±9% of galaxies in BvdHG have S24/S3 
color smaller than this value, while in the Karachentsev’s catalog, this ratio is only slightly 
lower, 89±14%, among the galaxies with same R band magnitude cutoff. That is, the small 
effect of the dispersion in infrared properties on the overall completeness is conhrmed by the 
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behavior of the complete sample of local galaxies. 

In Fig. 5, we plot the LFs obtained from the two spectroscopic catalogs along with the 
completeness limits. The hlled circles are resnlts from the BvdHC and the open squares from 
the MBC. The dotted vertical line is the detection limit at 24 /rm, the solid line is the limit 
of the MBC, and the dashed one is the limit of the BvdHC. The shapes of these two LFs are 
similar, but the LF from the MBC has an overall higher number density than that from the 
BvdHC. The reason for this difference is that Mobasher’s spectroscopic survey covered the 
whole cluster core, where the galaxy number density is the highest, but only a small portion 
of the outer region, while the BvdHC is based on a much larger area including both the core 
and the outer region. Also, the LF from the MBC has a larger variance compared to the LF 
from the BvdHC because it is based on a smaller sample. Below the 24 /xm detection limit, 
the loss of faint galaxies due to the limit of the IR observations causes quick drops in both 
LFs. Above this limit, the faint end slope of the MBC LF is steeper than the BvdHC LF, 
which is consistent with the different completeness limits of these two LFs. 

To have a more quantitative comparison, we fitted these two LFs with some analytical 
functions. We discarded all the data points below 10^^ ergs s“^ and use a chi-square mini¬ 
mization method to hnd the best htting parameters. Since we do not have many data points 
at the bright ends which are critical to determine it is important to constrain the htting 
beyond the last bin for the non-detection of brighter galaxies. In order to incorporate this 
factor into the htting, we calculate the integrated expected galaxy number brighter than the 
brightest galaxy actually observed for each trial function and use this number to estimate 
the probability of the non-detection. This integration is carried out from the brightest lu¬ 
minosity observed to a luminosity 2 orders of magnitude brighter. The results change little 
when we extend the integration to a higher upper limit. We include the chi-square of this 
non-detection into the total chi-square value for the minimization process. We hrst htted the 
LFs with the Schechter function (Schechter 1976) and the best-fitting parameters we found 
are: 


a = log(L;/Ls) = 10.48l»;«, tor the MBC; (1) 

a = 1.4VS:S1; log(i;/is) = 10.49^“;g, for the BvdHC. (2) 

The htting results are shown as the solid curves in Fig. 5. The L*^ values of these two LFs 
are very similar. The parameters of the MBC LF have larger uncertainties due to the few 
data points at the bright end to constrain the htting. It also has a steeper faint end slope 
than the LF for the BvdHC, which is in agreement with the fact that the MBC has been 
corrected for incompleteness at the faint end while the BvdHC has not. Considering this 
factor, we expect that the IR LF for the MBC gives a better estimate of the faint end slope 
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than the IR LF for the BvdHC. 

A recent work of Perez-Gonzalez et ah (2005) studies the 12 /xm LF from the 24 /xm 
emission of galaxies using the Spitzer data in two deep held surveys, the Chandra Deep Field 
South and the Hubble Deep Field North. Their results, coming from galaxies in the general 
held, provide a good comparison to our LF in a dense cluster. Their LF for galaxies with 0.0 < 
2 : < 0.2 gives the Schechter parameter of a = 1.23±0.07 and \og{L\ 2 /Lq) = 9.61T0.14. With 
a relation between Lir and L 12 given by Takeuchi et ah (2005), logLi^ = 1.02 + 0.972 logLi 2 , 
the obtained by Perez-Gonzalez et ah (2005) corresponds to a total IR luminosity of 
\og{L*^/L q) = 10.36+0.14, which is only slightly smaller than the value we got; the diherence 
is well within the one sigma error. However, we found a somewhat steeper slope at the 
low luminosity end (at about 2a signihcance). Perez-Gonzalez et al. (2005) suggest that 
incompleteness may have reduced the value of a in their LF. 

Rush, Malkan & Spinoglio (1993) obtained a LF using an all-sky 12 /im survey from the 
IRAS Faint Source Catalog, Version 2, and htted it with a two-power-law function 

t(L) = cn-"(i + ^)-'’. 

For the non-Seyfert subsample (the majority are normal galaxies, and about 5% are starburst 
galaxies), they found the best-htting parameters are a = 1.7, f3 = 3.6, and \og{L\ 2 lL q) = 
9.8. The average redshift of their non-Seyfert subsample is 0.013, comparable to our sample’s 
average redshift. Their result has a much steeper slope at the faint end compared to the 
Perez-Gonzalez et al. (2005) result and a little larger which corresponds to log(L*^/L 0 ) = 
10.55. To make a fair comparison to their result, we tried to £t our LFs with the same two- 
power-law function. However, since the two-power-law function has too many free parameters 
and our small sample size gives poor constraints at the high luminosity end, a free fitting 
failed to give a reasonable result. We hxed the slope index at the high luminosity end 
to the best-htting value given by Rush, Malkan & Spinoglio (1993) and kept a and 
the normalization free. The fitting gives very similar results to the Schechter fitting; the 
best-htting parameters are 

a = 1.481°;J2; log(L;/L 0 ) = 10.24l°;i; /3 = 3.6 {fixed), for MBC; (3) 

a = 1.38lg;}2; log(-^ir/-^o) = — 3.6 {fixed), for BvdHC. (4) 

Again, we note that the MBC LF gives a larger value for the faint end slope and the BvdHC 
LF gives a more reliable estimate for the L*. value. The MBC LF has a faint end slope a 
little hatter than Rush’s LF, and both LFs give L*. values smaller than Rush’s L*,) with 
signihcance about one sigma. 
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Takeuchi, Yoshikawa & Ishii (2003) estimated the 60 jum LF of the local galaxies in the 
Point Source Redshift survey of IRAS. Based on this LF, Le Floc’h et ah (2005) calculate 
the total IR LF using the 60 /rm total-IR IRAS correlation (e.g., Chary & Elbaz 2001) and 
£t the IR LF with a double-exponential function. The fitting gives a faint end slope index of 
1.23, log(L*^/L 0 ) = 9.25 and a = 0.72, where a is the parameter used to adjust the shape of 
the bright end of the LF (e.g., Le Floc’h et al. 2005). To compare to their result, we fitted 
our LFs with the same double-exponential function. Since neither of our LFs has enough 
data points at the bright end, we have to fix a = 0.72, the value they gave. By doing this, we 
found a faint end slope index of 1.52 ± 0.13 for the MBC LF and log(L*,/L 0 ) = 9.39 ± 0.37 
for the BvdHC LF. Again, the Coma IR LF has a very similar L*. to the held LF. The faint 
end slope of Takeuchi’s LF is similar to that of Perez-Gonzalez et al. (2005), and they both 
are flatter than the Coma LF. 

The comparisons between the Coma IR LFs and the IR LFs from the general held do not 
show significant variation of the L*^ value in different environments. However, a difference 
in the IR luminosity as small as the 0.3 mag ( 0.12 dex in luminosity) difference in of 
cluster LF and held LF shown by De Propris et al. (2003), is beyond the capability of our 
study. The faint end slope of the Coma IR LF is steeper than that of Perez-Gonzalez et al. 
(2005) and Takeuchi, Yoshikawa & Ishii (2003), but a little shallower than that of Rush, 
Malkan & Spinoglio (1993). This comparison, although complicated by completeness issues, 
does not support a strong dependence of the shape of the IR LF on environment. 

Despite the similarity in the shape of the Goma IR LF and the held IR LF, there might 
be a large portion of IR-inactive galaxies in the cluster compared with the held. Assuming a 
line-of-sight dimension of the Goma cluster of 13 Mpc, the BvdHG IR LF has a 0* value - the 
space density at - about 45 times larger than the 0* value given by Perez-Gonzalez et al. 
(2005). Because L*. is far above our detection limits and those of the spectroscopic surveys 
(using the proportionality we found for local galaxies between optical and IR luminosities), 
our study should be complete there. With the same assumptions, we hud that the average 
density at L\ is 62.9 ± 15.2 times that in the held (Geller et al. 1997; Beijersbergen et al. 
2002a). That is, the infrared-emitting galaxy density is only slightly less enhanced in the 
cluster than the optical galaxy density; there are few extra IR-inactive galaxies in the cluster. 
A second approach to this issue is to examine a sample of held galaxies, and see how many 
would be detected in the infrared using the same selection method as we have used in the 
Goma Gluster. To implement this approach, we have again used the catalog of galaxies 
within 10 Mpc from Karachentsev (2004). We have compiled the IRAS data for galaxies 
down to Mb = -17.5 mag, the completeness limit of the BvdHG. If these galaxies were at 
the distance of the Goma Gluster, we hud that 0.89 ± 0.12 would be detected above our 24 
pm limit, whereas the portion of infrared-detected galaxies in the BvdHG down to Mb = 
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-17.5 mag is 0.56 ± 0.05. In agreement with our first estimate, there is only a small dehcit 
of IR-active galaxies compared with the behavior in the held. 

It is possible that this difference is partly due to the morphology-environment correlation 
of galaxies because our Coma sample has about 53% early-type galaxies while Karachentsev’s 
held galaxy sample only has about 17% early-type galaxies. So we divided our sample into 
early- and late-type subsamples and calculated the portion respectively. It turns out this 
portion just slightly increases in the late-type subsample compared with early-type subsample 
(57% vs. 55%). Therefore, the morphology-environment correlation can not account for the 
diherent portion of IR-active galaxies in Coma and the held. 


6.4. Luminosity Function in Different Regions of the Cluster 

Given the large coverage of the BvdHC, we can study the LFs in different regions of 
the cluster. Although the BvdHC is only complete down to R = 16.27, our previous results 
show the incompleteness probably will only make the faint end slope a little shallower. 
Since the completeness does not change very much across the cluster, incompleteness will 
not bias the comparison of LFs in the different regions. Following Beijersbergen et al. 
(2002a), we define the core region of Coma as the area with r < 12' (~ 0.3 Mpc) centered 
at a = 12^59™'43®, S = -|-27°58'14". We also define an annulus region outside of the core 
for 12' < r < 24' (~ 0.6 Mpc). Another interesting area is the group of galaxies around 
NGC 4839. It is the second densest region in the cluster and its X-ray emission suggests 
that the group is falling into the cluster (Neumann et al. 2001). The interaction between 
the group and the cluster may trigger star forming activities and therefore affect the IR LF. 
Therefore, we also select the circular region centered at NGC 4839 with the same radius as 
the Coma core region. We constructed LFs in these three regions. Apart from these regions, 
we took the rest of the area with MIPS coverage as a whole to be the outskirt region of 
the cluster. A sky map of these regions and all the galaxy members detected at 24 /rm is 
shown in Fig. 1. There are 40, 56 and 28 galaxies detected in the IR in the Coma core, 
the surrounding annulus and the NGC 4839 region. In the outskirt region, 101 galaxies are 
detected. The two circular regions have an area of 0.13 deg^, the annulus region has an area 
of 0.38 deg^, and the outskirt has an area of about 2.41 deg^. Therefore, the ratios of the 
projected number density of infrared emitting galaxies in the Coma core, NGC 4839 region, 
annulus, and outskirt region are about 6:4:3:1. Even the lowest-density region has a space 
density of infrared galaxies roughly 40 times that in the held. 

All the LFs were fitted with the Schechter function. The results are shown in Fig. 6 
and the best fitting parameters are listed in Table 1. Because the small number statistics in 
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small regions cause large uncertainties in the Schechter function fitting, simple comparisons 
of the best htting parameters between these LFs are ambiguous and need to be taken with 
care. The large uncertainties at the bright end of the LFs may cause very different htting 
results for the exponential cut-off of the Schechter function and, therefore, unreliable L* 
values. However, interesting variations are apparent in other aspects of the LFs. The Coma 
core region has a hatter faint end and fewer luminous galaxies compared to the LFs in the 
other regions. In particular, when we compare it with the LF in the annulus region, which 
has similar total number counts as in the core, it is apparent that galaxies in the core region 
are lacking at the high luminosity end. A hatter faint end in the core indicates a lack of faint 
galaxies as well. The NGC 4839 and the annulus regions have similar number densities and 
their LFs are not very diherent from each other at the faint end. They both have a steeper 
faint end LF than the Coma core. At the bright end, it seems that the annulus region has 
more luminous galaxies than the NGC 4839 region does, but the diherence is not signihcant 
given the uncertainties. The LF in the outskirt region is better constrained at both faint 
and bright ends due to its larger number of galaxies. Its faint end slope is steeper than 
the Coma core region but shallower than the annulus and NGC 4839 region. However, this 
faint end slope is largely constrained by the lowest point in our htting process (the point at 
Lir = ergs s“^) and we suspect that incompleteness may have a more severe effect on 

this point in the outskirt region than in the other regions (e.g., because of the lower density 
of cluster members on the sky). If we discard this point in our htting, we have a much 
steeper faint end slope with a = 1.52^q;} 7. With this correction, there appears to be a trend 
of steeper faint end slope toward the outer regions of the cluster, similar to the behavior in 
the optical bands (Beijersbergen et al. 2002a; Beijersbergen & van der Hulst 2003). Although 
the faint end slope may be questionable, the bright end of the LF in outskirt region is well 
constrained and has a very similar L* value as the total LF. In summary, we found that the 
Coma core region lacks both very faint and very bright galaxies compared with the outer 
regions. The NGC 4839 region does not show signihcant difference in the LF from that of 
the annulus region with similar number density. It is also worth noting that all the galaxies 
with Lir > 10^"^ ergs s“^ reside outside of the core region. 

All the LFs in the Coma cluster calculated above are actually the 2-D projection of the 
real LF, and therefore will be affected by projection effects (Valotto, Moore & Lambas 2001; 
Beijersbergen et al. 2002b). The projection effect is most serious in the core region and it will 
probably make the faint end slope steeper. That is, the hat faint end slope of the core region 
may become even hatter after deprojection. On the other hand, the projection ehect will 
make the lack of bright galaxies in the core region more severe. Fig. 1 also shows the galaxies 
with Lir > 10^^'^ ergs s“^ as star signs. They are more or less uniformly distributed in the 
whole region, without any concentration in the core region or the NGC 4839 region. The 
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deprojection from this 2-D distribution will make the bright galaxies move further outward. 

The interpretation of the galaxy population variations across the cluster only from the 
change of the shape of the LFs may be misleading without knowledge of the fraction of 
galaxy members detected in the IR. We already know that the overall fraction of optical 
galaxy members detected at 24 /rm is about 44% for the BvdHC. This fraction is smallest 
in the core region, at about 37 ± 7%, and it is about 46 ± 4% in the outer region. The 
difference is not very significant considering the large statistical errors. If we only consider 
the fraction for the galaxies brighter than the completeness of the BvdHC, the difference 
is even smaller, with the fraction about 54 ± 9% in the core and 57 ± 5% in the outer 
region. Also, a correction for the projection effect, if possible, would make the difference 
even smaller. Thus, the fraction of the galaxy members detected in the IR does not change 
very much across the cluster, providing a uniform foundation for the comparison of the shape 
of the IR LFs. 


6.5. Contribution of the Different Types of Galaxies to the Total LF 

Since the BvdHC also has morphology information for each galaxy, we can study the 
contribution of early type (E/SO) and late type galaxies to the total LF. In Fig. 7, we plot the 
LF of the late type galaxies and early type galaxies along with the total LF for comparison. 
The late type galaxies here include all spirals and irregulars. 

From Fig. 7, we find that the early type galaxies make a larger contribution to the 
number counts of the LF than the late type galaxies at < 10"^^ ergs s“^, while the 
late type galaxies dominate the bright end of the total LF. This behavior indicates that 
although the late-type galaxies dominate the bright population, there are more faint early- 
type galaxies than faint late-type galaxies. However, we note that at the faint end, it is 
possible to misidentify a spiral as an SO galaxy and, therefore, the number of early-type 
galaxies may be overestimated. We ht the LF of late type galaxies with the Schechter 
function, and the best-htting parameters are 

a = i-isiS:!; log(i;/ia) = 10.44«:2 (5) 

This LF has a flat faint end and a similar L*^ to the total LF. Although the a we derive 
would have a higher value if the incompleteness were taken into account, it is still smaller 
than the index of the total LF that is affected by the incompleteness in the same way. The 
steeply rising faint end of the total LF is boosted by the increasing number of early type 
galaxies with low IR luminosity. 

Using ISO data, Pozzi et al. (2004) deduced the 15 /rm LF of the European Large Area 
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ISO survey (ELAIS). The index of the faint end slope they found for the spiral galaxies with 
2 : < 0.2 is very close to our value, with a = 1.10 ± 0.25. 

Using the LFs of early type and late type galaxies, we can calculate the surface density 
of the total IR luminosity of these two groups down to the detection limit of the 24 /rm 
observations. It turns out that the surface density of IR luminosity contributed by early 
type galaxies is only about 15% of the total surface density. Therefore, the early type 
galaxies make a rather small contribution to the total IR luminosity of the cluster, but they 
make a signihcant contribution to the number counts of faint galaxies and therefore affect 
the shape of the LF. 


6.6. Measuring the SFR from the IR LF 

Since IR luminosity is a good tracer for star-forming activity, the IR LF allows us to 
estimate the total SFR of the cluster. Although this paper is the first work reaching such a 
depth in the IR luminosity of the Coma galaxies, which means detecting a lower level of star 
forming activities, there are other works measuring the SFR at a higher level by measuring 
the ionization lines. Iglesias-Paramo et al. (2002) used a deep wide held Ho survey of the 
Coma cluster to deduce the Ha LF. They detected 22 sources in the Ha band. Five of them 
are not in either the BvdHC or the MBC. We detected all the rest of them at 24 fim and 
hence obtained the for them. Using the conversion formula given by Kennicutt (1998b) 

SFR{MqYt-^) = 4.5 X 10-^^Li^(ergs s"^) (6) 

we deduce the SFR for these objects and compare them with the SFR given by Iglesias- 
Paramo et al. (2002). The result is shown in Fig. 8 . The two results are basically consistent, 
but the SFRs measured from Lir are larger on average, and this discrepancy is more pro¬ 
nounced in the galaxies with higher SFR. This discrepancy was also found by Kennicutt 
(1998a) when he compared the SFR deduced from the measurement of Lir and the Bry 
emission line. He justihed the discrepancy by citing effects of extinction and the heating of 
dust by longer lived stars than those exciting the emission lines. These arguments are also 
applicable to our case. 

The Ha survey has a smaller area coverage than the MIPS observations and is much less 
complete than the IR survey. For the 18 objects detected with both Ha and 24 /rm emission, 
16 have Lir > 10^^ ergs s“^, but in our sample we have 40 sources with Lir > 10^^ ergs s“^. 
Using the Ha LF, Iglesias-Paramo et al. (2002) showed the SFR density of the Coma cluster 
to be 1.36 M 0 yr“^Mpc“^ by integrating the best htting function over the whole range of 
luminosities and assuming the radius of the Coma cluster to be 6.5 Mpc. In our case, we 
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calculate the total IR luminosities from the best htting Schechter functions in the range of 
1042 gj^gg g-i ^ ^ ]^g 45 gj,gg g-i_ rjj^g uppgj; limit corresponds to the brightest galaxies 

actually observed and the lower limit excludes the part of the LF with serious incompleteness. 
The formula we used to convert the IR luminosity to SFR will also become problematic below 
this lower limit. Fortunately, from the shape of the LF, the galaxies in the 10^^ ergs s“^ < 
Lir < 10 ^^ ergs s“^ luminosity range dominate the total luminosity and there is only a small 
difference between our integration and an integration over the whole luminosity range. The 
total SFRs of the cluster, deduced from the best htted Schechter functions of the MBC and 
the BvdHC IR LFs, are 24.48 M 0 yr“^Mpc“^ and 11.41 M 0 yr“^Mpc“^ respectively. The 
lower value given by the BvdHC LF is because the BvdHC covers more outskirt regions 
where the IR galaxy densities are relatively smaller. The BvdHC is also less complete than 
the MBC, and thus underestimates the contribution from the faint IR galaxies. The total 
SFR in the 8.5 Mpc^ area of the BvdHC survey and MIPS observation is about 97.0 

If we assume the region we observed in the IR has the same line-of-sight dimension 
as assumed by Iglesias-Paramo et al. (2002), e.g., ~ 13 Mpc, these two IR LFs give SFR 
densities of about 1.88 and 0.88 . The SFR density deduced 

from the Ha LF is smaller than the value given by MBC but larger than that given by 
BvdHC. However, considering the Iglesias-Paramo et al. (2002) survey also covers mostly 
the central region, we calculate the SFR from the BvdHC in the region of the Ha survey 
coverage and obtain a larger SFR of 1.7 M 0 yr“^Mpc“^. Therefore, the IR LF gives a more 
complete estimate of the total SFR of the cluster than the available Ha LF. If we adopt 
0.88 M 0 yr“^Mpc“^ as the general SFR density of the Coma cluster, we hnd it to be about 
60 times larger than the SFR density of the general field (Perez-Gonzalez et al. 2005). This 
difference is comparable to the difference in (p* between these two LFs. Thus, the higher 
SFR density in the Coma cluster is mainly due to the overall higher IR galaxy density in 
the cluster than in the held, not to any differences in the shape of the IR LF. 


7. Discussion 

Despite the evidence that the SFRs of galaxies are different in cluster and held regions 
(Gomez et al. 2003; Balogh et al. 1998), the comparison of our IR LFs of the Coma cluster 
with the IR LFs of held galaxies does not support a strong dependence of the shape of the 
IR LF on environment. However, the measurements of the SFR in Gomez et al. (2003) 
and Balogh et al. (1998) are based on the bright galaxies {R < 17 mag) with ionizing 
emission lines, which approximately corresponds to the galaxies with L^r > 10 *^^ ergs s“^. 
In fact, Gomez et al. (2003) found that the correlation of the SFR and the environment is 
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most noticeable for the strongly star-forming galaxies. So, it is very possible that our IR 
LF, because it is not well constrained at the bright end, does not show a difference from 
the held LF simply because of the lack of enough very luminous IR galaxies in the small 
sample to draw any meaningful conclusions. In any case, the similar of the Coma cluster 
to that of the general held LF, as well as the fairly good ht with similar functions up to 
Lir ~ 10^^ ergs s“^, are evidence against a strong correlation between the global SFR for 
infrared-bright galaxies and their environment. 

Explanations for this behavior may come from some recent works by Balogh et al. 
(2004a,b). Balogh et al. (2004a) studied SFRs of galaxies in group and low-density envi¬ 
ronments and found that although the fraction of the star-forming galaxies is very sensitive 
to the galaxy density, the distribution of IFo(Hq,), the equivalent width of Ho;, in the star¬ 
forming galaxies is independent of environment. Balogh et al. (2004b) studied the color 
distributions of galaxies in diherent environments including the typical environment of a 
cluster core. After dividing their sample into red and blue components in several luminosity 
bins, they found that the ratio of these two components is a strong function of the local 
density but the mean value and the shape of the color distribution of each component are 
nearly independent of environment. They proposed that most star-forming galaxies today 
evolve independently from their environment. Both these works suggest that interactions of 
the galaxies, probably happening in a very short time scale (r < 0.5 Gyr), may be respon¬ 
sible for triggering star formation and leave the galaxies to evolve afterwards independently 
from their environments. Our results appear to agree well with these works and push the 
independence of the global SFR of IR galaxies on environment to an even lower level of SFR. 

It is also possible that the lack of difference between the shape of the Coma LF and 
LFs from the general held results in part from the effects of averaging. First, the general 
held also includes galaxies in clusters. Secondly, and more importantly, the galaxies in the 
outskirt region have evolved largely in the held and are only beginning to fall into the cluster. 
However, their LF of dominates the total LF of the cluster at luminosities above L*^. Since 
the outskirt LF has a similar L*. and a little smaller a compared with the total LF, it is 
possible that the contribution from the dense region mostly goes to the faint end of the LF 
where we do detect an environmental ehect. The early-type galaxies also contribute greatly 
to the number counts in this range. However, exact comparisons are complicated by the 
incompleteness. It has been reported that there is a group of dwarf galaxies contributing 
to the steeply rising faint end slope of the optical and near-IR LF in the Coma cluster 
(Bernstein et al. 1995; Seeker & Harris 1996; Trentham 1998; Mobasher & Trenthan 1998; 
De Propris et al. 1998). These dwarf galaxies are beyond our completeness limit {R < 19 
mag) and are not the contributors of the steeply rising faint end of the IR LF. 
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Although we did not obtain any direct evidence of a correlation of SFR and environment 
in the comparison of the shape of the Coma IR LF and the general field IR LF, there is some 
evidence of a change in SFR with the environment inside the cluster. We found that although 
the fractions of the galaxy members detected at 24 /rm do not change very much across the 
cluster, the LF in the core region has a flatter faint end and is lacking the bright IR galaxies 
compared to the LF in the outer region. In addition, all the galaxies with L*. > 10"^^ ergs s“^ 
lie outside of the core region. This behavior shows that the strongest star-forming activity 
happens in the lower density region of the cluster. This is consistent with the speculation 
that galaxies in a crowded environment lack gas and dust due to the interactions between 
the galaxies or between the galaxies and the cluster potential well and therefore can not 
support a large SFR. The flatter faint end in the core region, on the other hand, suggests a 
deficiency of faint IR galaxies as well. 

Mori & Burkert (2000) studied the gas stripping of dwarf galaxies by ram pressure in 
the cluster and found that dwarf galaxies will lose virtually all of their gas instantaneously 
if their core mass is smaller than a critical mass {Mcr)- For the Coma cluster, \og{Mcr/M q) 
is about 10.9 at the median distance of the cluster galaxies. This value increases to about 
11.7 and 12.9 at r = 24' (~ 0.6 Mpc) and r = 12' (~ 0.3 Mpc). Mori & Burkert (2000) 
would argue that galaxies with core masses smaller than these values will lose all of their gas 
very quickly (r ~ 10® yr). Thus, even if the triggering of the star-forming activities in these 
galaxies happens in the galaxy group before they fall into the cluster core, gas stripping will 
prohibit them from keeping up such activities by depriving them of fuel. 

However, such an effective stripping contradicts our result. The faint IR galaxies in 
our LFs (41.9 < \ogLir < 42.6, about the range of the lowest two points in the LFs above 
the incompleteness) have R ^ lb mag on average. The early type galaxies in the Coma 
cluster usually have a stellar mass-to-light ratio smaller than 8 (Jprgensen 1999). The stellar 
mass is comparable to the core mass of a galaxy, so we can use this ratio to estimate the 
core mass of a galaxy from its luminosity. Considering that faint galaxies usually have 
relatively large M/L, we take M/L ^ 20 as a conservative upper limit for these galaxies. We 
also assume all cluster galaxies have similar M/L ratios. With this ratio and the average 
R magnitude of the faint IR galaxies, we obtained the upper limit of the mass for these 
galaxies as \og{M/M q) k, 11.0. This mass is smaller than the critical mass for gas stripping 
at r < 24' (~ 0.6 Mpc), which means most of the faint galaxies in this region should already 
have lost all of their gas, so we would expect a drop in the number counts of the IR LF. 
However, the IR LF in the annulus region does not show such a drop, although we do see a 
flattening of the slope in the core region. If this flattening of the faint end slope is the result 
of the total gas stripping, there seems to be a factor of ten discrepancy between the critical 
mass Mori & Burkert deduced and the one suggested by the IR LFs. A smaller critical mass 
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for these IR galaxies is also consistent with the small difference we fonnd in the Coma IR 
LF and the IR LF of general held. Mori & Bnrkert (2000) also pointed out that there were 
some issues they did not consider in their simulation that may affect the critical mass, e.g., 
the heating of the gas from star forming activities. For the IR galaxies, this effect may be 
very important because it would lower the critical mass and provide a possible explanation 
of why these faint galaxies still have star forming activities. An underestimated M/L ratio 
can also contribute to this difference. However, to resolve the discrepancy in this way, the 
M/L would need to be as high as 100, which seems unlikely. 

We also observed a high ratio of early-type galaxies in the core region of the Coma 
cluster, ~ 80%; this ratio drops to about 50% for the other regions. The change of the shape 
of IR LF in the core region is also possibly caused by a morphology-environment correlation 
rather than a SFR-environment correlation. 


8. Conclusions 

Using MIPS 24 /rm observations and two spectroscopic surveys of the Coma cluster, 
we present the IR LF of the cluster. The shape of the Coma cluster LF does not differ 
from that of the general held signihcantly. The L*. value of our LF is very similar to those 
given by Rush, Malkan & Spinoglio (1993) and Perez-Gonzalez et al. (2005), which are both 
based on surveys of general helds. The faint end slope of the Coma cluster is shallower than 
the slope of Rush, Malkan & Spinoglio (1993) but steeper than that of Perez-Gonzalez et al. 
(2005) and Takeuchi, Yoshikawa & Ishii (2003), again indicating little variation between held 
and cluster. In addition, the overall proportion of IR-active galaxies in the cluster is only 
slightly less than in the held. Thus, the overall pattern of star formation in cluster members 
is surprisingly similar to that in the held, despite an increased galaxy space density by an 
average factor of 40. 

However, in the cluster core where the galaxy density is six times higher still, we found 
a shallower faint end slope and a smaller L*. compared to the outer region of the cluster, 
which indicates a decrease in the number of faint IR galaxies as well as in the very bright 
ones. The IR-bright galaxies are distributed around the outer region of the cluster. All the 
galaxies with > 10^^ ergs s“^ lie outside of the core region, e.g, r > 340 kpc. No special 
feature of the IR LF was found in the NGC 4839 region. 

In determining the LF of diherent morphological types, we found that early type galaxies 
only make about a 15% contribution to the total IR luminosity density, but they dominate 
the number density at the low luminosity end. The global SFR density in the cluster is 
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about 0.88 MqJt ^Mpc ^ and the total SFR in the 8.5 Mpc^ area of the central cluster is 
about 97.0 
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Fig. 1.— The sky map of the Coma cluster. The dots are the galaxy members detected 
at 24 pm in the BvdHC. The stars designate the sources with ergs s“^. The 

circle in the center dehnes the Coma core region, with r < 0.°2, and the surrounding annulus 
is dehned as 0.°2 < r < 0.°4. The circular region southwest of the Coma core is centered 
at NGC 4839 with same radius as the Coma core. The dotted rectangular region is the 
MIPS 24 pm coverage. The dashed region is the survey region of the BvdHC, and the two 
dash-dotted rectangular regions are of the MBC. 
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Fig. 2.— The color - color/flux plot of the early- and late- type galaxies. The open circles 
and open triangles represent early- and late- type galaxies, respectively, (a) the 70 - 24 color 
vs. 24 - R color. The crosses are the color of the template SEDs. The dotted line is the 
detection limit set by the completeness of the 70 pm observations, (b) the 24 pm flux density 
vs. 24 - R color. The galaxies also detected at 70 pm are plotted as hlled symbols. 
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Fig. 3.— The Li^jL 2 A, ratio as a function of 24 - R color. The diamonds denote the normal 
spirals, the square a LIRG and the triangles designate ULIRGs. The dots are the results of 
the interpolation of the Goma galaxies from their 24 - R colors. The Li^jL 2 i ratio given by 
(Lagache et ah 2003) for a normal spiral is shown as the dashed line. 
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Fig. 4.— Comparisons of obtained from different color correlations. The filled circles are 
the galaxies with S '70 greater than the 80 mJy completeness limit, and the open circles are 
the galaxies under the limit. 
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Fig. 5.— (a) The IR luminosity function of the Coma galaxies. The hlled circles and the 
open squares represent the LF from the BvdHC and the MBC, respectively. The error bars 
denote the statistical error. The three vertical lines are completeness limits: the dotted 
vertical line shows the completeness of the 24 /rm detections, the solid vertical line shows 
the completeness of the MBC after correcting so far as possible for incompleteness, and the 
dashed vertical line shows the completeness of the BvdHC. The solid curves are the results 
of htting the LF with the Schechter function and the dashed curves are the results of htting 
with a two-power-law function with a hxed slope at the bright end. (b),(c),(d) and (e) are 
the error contours for the htting parameters L*^ and a. The contour levels are 1, 2 and 3 a. 
The best-htting parameters are indicated by the cross symbols. 
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Fig. 6.— The LFs of the galaxies from the BvdHC in different regions of the Coma cluster 
and the corresponding 1, 2, 3 cr error contour maps for the Schechter function fitting. The 
vertical lines are the same as in Fig. 5 and the solid curves are the best fitting Schechter 
functions. The best-htting parameters are indicated by the cross symbols on the error contour 
maps. 
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Fig. 7.— The contributions of the early- and late- type galaxies to the LF. The galaxies are 
all from the BvdHC. The vertical lines are the same as in Fig. 5. The total LF and the LF 
of the late-type galaxies are htted with the Schechter function, and the results are shown as 
the solid curves. 




Table 1. 


The best fitting parameters of the LFs in different regions. 


Region 

log(-^ir/ -^ 0 ) 

a 

NGC 4839 r < 12' 

1 n 21 

iU.Zi_o 48 

1 44+0-30 

core r < 12' 

Q 44+0-27 

n qq+0.34 
'O-^y-0.52 

12' < r < 24' 

lo.eslJI? 

1 42 + 0-20 
-‘-•+^-0.20 

out skirt 

10.50l°j° 

1 Qo+ 0.11 
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